Spin States in Graphene Quantum Dots 
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We investigate ground and excited state transport through small (d ~ 70 nm) graphene quantum 
dots. The successive spin filling of orbital states is detected by measuring the ground state energy as 
a function of a magnetic field. For a magnetic field in-plane of the quantum dot the Zemann splitting 
of spin states is measured. The results are compatible with a g-f actor of 2 and we detect a spin- filling 
sequence for a series of states which is reasonable given the strength of exchange interaction effects 
expected for graphene. 
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Spin qubits in quantum dots pQ are interesting can- 
didates for the implementation of future quantum infor- 
mation processing. Single spin preparation, manipula- 
tion and read-out has so far been demonstrated predom- 
inately in GaAs-based systems [2j [3]. Spin coherence 
times in such systems are limited by hyperfine coupling 
and spin-orbit interactions. In graphene-based nanos- 
tructures both limitations are expected to be significantly 
reduced in strength j2HZ]. The electostatic tunability of 
graphene quantum dots as well as the observation of ex- 
cited states have been demonstrated recently [5UT2] . The 
g-factor in graphene has been measured via conductance 
fluctuations [13] . but the general spin properties of con- 
fined electrons in graphene have remained elusive. 

Here we present transport measurements on a 
graphene quantum dot as a function of magnetic field 
and analyze the evolution of Coulomb peaks and excited 
states in both perpendicular and parallel magnetic fields. 

A scanning force microsope (SFM) picture of the sam- 
ple is shown as an inset in Fig. 1(a). The graphene quan- 
tum dot in the center is connected via tunneling constric- 
tions to source (s) and drain (d) contacts. It is tunable by 
the back gate (bg) and several in-plane gate electrodes, 
including the plunger gate (pg). The fabrication of the 
sample is described in detail in Ref. [T4J All measure- 
ments were performed in a dilution refrigerator with a 
base temperature of 90 mK. 

The source-drain conductance of the device measured 
over a large back gate voltage range is shown in Fig. 1(a). 
The region of suppressed conductance between and 
13 V ("transport gap" [15]) separates the hole-transport 
regime at low Vbg from the electron transport regime 
at high Vbg. Coulomb-blockaded transport through the 
quantum dot is observed in this region at closer inspec- 
tion. Regions A and B indicate the gate voltage ranges of 
the two measurement regimes investigated in this paper. 
Previous measurements in region A [16] gave evidence 
for an electron-hole crossover in the Coulomb-blockaded 
quantum dot around Vb g = 5.5 V. 

The same conductance resonances scanned with the 
back gate in region A (B) can also be investigated by 
keeping Vbg fixed at -0.9 V (-1.8 V) and sweeping the 
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FIG. 1: (Color online) (a) Quantum dot source-drain con- 
ductance for varying back gate voltage (Vt> = 100 /iV). The 
inset shows a scanning force microscope picture of the mea- 
sured graphene quantum dot device, (b) Coulomb diamond 
measurement recorded on the hole side of the electron-hole 
crossover. The gate dependence at Vb = 11 mV (white dashed 
line) for increasing perpendicular magnetic field is plotted in 
(c). Both edges of the diamond in (b) show similar magnetic 
field dependences typical for spin pairs. 



plunger gate voltage [see top-axis in Fig. 1(a) for the cal- 
culated plunger gate voltage at V\> g = —0.9 V (region A) 
corresponding to the measured back gate voltage. This 
plunger gate scale has been used throughout the paper 
as a gate voltage scale to ease comparison] . 

Two Coulomb-blockade resonances from regime A are 
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investigated by measuring the differential conductance as 
a function of plunger gate voltage and source-drain bias 
at B = T in Fig. 1(b). At V h = mV the differential 
conductance is suppressed except for the Coulomb reso- 
nance conditions (circles) where two conductance peaks 
A2 and A3 are observed. By increasing the bias window 
a diamond shaped region of blocked conductance is mea- 
sured, revealing an electron addition energy of 20 meV. 
The change in differential conductance parallel to the 
edges of the diamond is attributed to excited states pro- 
viding additional transport channels. Their energy split- 
ting from the ground state is obtained directly from the 
measurement (e.g. 6 meV for the state marked with the 
gray arrow). 

A cut through the diamond at fixed bias voltage = 
11 mV [dashed white line in Fig. 1(b)] for increasing 
magnetic field perpendicular to the graphene plane is 
shown in Fig. 1(c). After a kink around B = 3.5 T the 
edges of the diamonds move towards higher gate voltages 
with increasing field. This kink has been interpreted in 
Refs.[16j[T7las an indication of filling factor two for holes 
in the quantum dot. The same kink is observed for the 
excited state resonance at V pg = 3.66 V. The comparable 
evolution of the two diamond edges (white arrows) could 
be a result of spin pairing, i.e., two subsequently filled 
electrons occupying the same orbital state with opposite 
spin orientation. 

Such pairs are interesting candidates for studying spin 
properties. In general, the extraction of a small Zeeman 
splitting from the magnetic field dispersion of individ- 
ual conductance resonances is difficult, because any small 
orientational misalignment of the sample causes a signifi- 
cant perpendicular magnetic field component. Typically, 
the dispersion of states in the perpendicular field domi- 
nated by orbital effects is much stronger than the faint 
Zeeman splitting. Hence, small misalignment severely 
hampers the extraction of precise g-factor values from 
the dispersion of individual peaks. The problem can be 
circumvented by analyzing the peak-to-peak spacing of 
spin-pairs, because orbital contributions can be signifi- 
cantly reduced by subtracting the positions of individual 
peaks sharing the same orbital shift p~8j [19] . 

Potential spin pairs are also identified in the gate 
regime B. In Fig. 2(a) the evolution of Coulomb peaks 
with increasing perpendicular field is displayed. The low- 
est two peaks (Bl, B2) and the following two (B3, B4) 
are identified as potential spin pairs due to their similar 
peak evolution. In Fig. 2(b) a measurement of the same 
peaks is shown after careful in-situ rotation of the sample 
into an orientation parallel to the magnetic field. 

To analyze the movement of the peaks shown in Fig. 2, 
the peak positions are extracted by fitting. The results 
for the two pairs (B1-B4) in Fig. 2 and the pair (A2-A3) 
in Fig. l(b,c) are shown in Fig. 3(a). Peak positions are 
plotted with suitable offsets in V pg such that pairs coin- 
cide at B = T. For low perpendicular magnetic field 
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FIG. 2: (Color online) (a) Coulomb peaks as function of per- 
pendicular magnetic field recorded at V& = 100 fiY. The con- 
ductance is plotted in linear scale to enhance visibility of 
amplitude changes, (b) Same Coulomb peaks measured in 
parallel (in-plane) magnetic field after tilting the sample by 
90°. Here the conductance is plotted in logarithmic scale to 
increase the visibility of the peaks. 



(B± < 3 T for A2/A3, B1/B2; B± < 6 T for B3/B4) 
the orbital states of each pair have approximately the 
same magnetic field dependence. Hence spurious orbital 
contributions to the peak spacing resulting from slight 
misalignment are limited. The data in Fig. 3(a) for par- 
allel orientation of the magnetic field shows roughly linear 
splitting of pairs with increasing B\\. We interpret this 
splitting as the Zeeman splitting of two spin states with 
identical orbital wave functions. 

The magnitudes of the Zeeman splitting of the three 
pairs discussed above and for two additional peak- 
spacings from region A are plotted in Fig. 3(b). The 
data is obtained in two steps. First, the peak-to-peak 
spacing from Fig. 2(b) and corresponding measurement 
in regime A is extracted. Second, the peak spacing is lin- 
early fit from 4.2 T to 13 T and the offset obtained from 
the fit is subtracted. The Zeeman spin splitting AE Z is 
given bv [T91I20] 

A4+i = (S N +i - 2S N + S N -i) g^B\\ + const. 

Here Sn is the spin quantum number along the B\\- 
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FIG. 3: (Color online) (a) Comparing the evolution of three peak pairs in perpendicular (left) and parallel (right) magnetic 
field, (b) Coulomb peak spacing as function of B\\ for the three pairs shown in (a) and the two subsequent peak spacings 
(A4-A3 and A5-A4) in regime A. The dashed lines show the Zeeman splitting AE Z = ±\g\[iBB for a g- factor \g\ = 2. All 
peaks are converted to energy using the lever arm a pg = 0.13 ± 0.01 as extracted from the corresponding Coulomb diamond 
measurements, (c) Slopes obtained by a linear fit for B\\ — 4.2 — 13 T of peak-to-peak spacings shown in (b) as a function of 
Vp g . The slopes of measurement B are compensated for a parallel magnetic field misalignment of /3b = 3° (see text), where the 
shift induced by the compensation is depicted by the gray arrows. The peak difference B7-B6 has a strong kink at B\\ ~ 9 T 
and is fitted only up to 9 T. The spin- filling-sequence on top is extracted from the points along the dotted gray line. 



direction of the dot with N electrons. Spin differences 
of successive ground states A5^ x = Sjv+i — Sn are 
half-integer values (e.g. A^ 1 ) = 1/2 for adding a spin 
up electron or AS^ 1 ) = 3/2 for adding a spin up electron 
while "flipping" another spin from down to up). Hence 
differences between three successive spin ground states 

(2) 

take on integer values ASj^ +1 — Sjv+i — 2Sn + Sn-i — 
0, ±1, ±2, .... Apart from the splitting of the pair B2- 
Bl the splitting of all pairs in Fig. 3(b) show slopes 
AE Z = ASWgti B B\\ with integer values AS^ = 0, ±1 
and a g-factor value of approximately 2 [black dashed 
lines in Fig. 3(b)]. 

We now attempt to extract spin-related shifts from in- 
dividual unpaired conductance resonances despite the in- 
herent difficulty mentioned above. The final goal is to re- 
construct a plausible spin-filling sequence of the quantum 
dot in regimes A and B. We apply the following strategy: 
first, the misalignment angle j3 is estimated, such that 
second, the peak shifts in parallel field can be corrected 
for the orbital effects arising from the misalignment. Dif- 
ferent angles (3 exist for measurements in regimes A and 
B, because the sample had been rotated into perpendicu- 
lar orientation between the measurements B and A. The 
misalignment angle j3 was estimated assuming g = 2 and 
linear peak shifts in B± in the range < B± < 2 T and 
B\\ in the range 4.2 T < B\\ < 13 T. The angle /3 was then 

chosen to generate approximately AS^^ = ±1/2 for all 
peaks considered. The obtained values are Pa = 0°±0.5° 
and /3b = 3° ± 1°. As a consequence no correction for 
orbital effects was necessary in regime A. 

The misalignment-compensated-slopes of the peak- 



spacings are shown in Fig. 3(c) as a function of the 
plunger gate voltage V pg . The gray arrows indicate the 
change of the slope by the compensation of the angu- 
lar misalignment. The error bar includes the error from 
both the lever arm estimation (10% of the slope) and 
the variance from the linear fit of the peak-to-peak spac- 
ings. Three dashed black lines indicate the slopes cor- 
responding to AE Z /AB = {0, ±1} x 2/iB and the data 
points are colored according to their closest A^^ 2 ^ value 
(blue for AS< 2 ) = -1, green for AS< 2 ) = and red 
for AS't 2 ) = +1). Note that the two regimes A and B 
overlap: pairs B2-B1 and A5-A4 are in fact the same 
Coulomb peaks as it can be seen from a conductance 
measurement in dependence of V pg and Vb g (not shown) . 
We can extract the spin-filling-sequence ittiltti from 
peak A2 to B6 (dotted gray line). We emphasize that 
we see a clear deviation from titl as it has been seen 
in the low carrier regime of carbon nanotube quantum 

dots ethh. 

This result can be made plausible by considering the 
exchange interaction between the charge carriers. The 
strength of the Coulomb interaction in graphene can be 
characterized by the ratio r s between interaction energy 
and kinetic energy [25] 

Eint e 2 /(47re e r r ee ) 

r s = = — = a g ~ 0.9, 

^/7rh 2 v^n s 

independent of the charge carrier density. Here n s is 
the density of states, r ee = (7rn s ) _1//2 is the mean elec- 
tron separation, the Fermi velocity in graphene and 
a g = a ■ vf/€ y c the graphene fine- structure constant. For 
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FIG. 4: (Color online) (a) Close up of a Coulomb diamond 
measurement around the electron- hole crossover at B\\ — T. 
(b,c) Cross section at Vb — 12 mV [white line in Fig. 3(a)] 
for increasing in-plane magnetic field. In (b) the traces are 
offset in conductivity and smoothed over a five point window 
for clarity. The dashed lines in (c) have two different slopes 
corresponding to a Zeeman spin splitting with g — 2 and are 
plotted for comparison over the raw data. 



Zeeman energy AE Z = \g\n^B with g = 2 are plotted on 
top. Despite a slight underestimation of the splitting, the 
lines agree reasonably well with the measurement. Hence 
we attribute the two slopes to Zeeman spin splitting. The 
absence of clear peak splittings that would be induced by 
lifting the spin degeneracy with B can be seen as another 
indication of significant charge carrier interactions. 

In summary, we have analyzed a graphene quantum 
dot in an in-plane magnetic field by investigating the con- 
ductance peak spacing and the evolution of excited states 
around the electron-hole crossover. In both regimes A 
and B spin-pairs were identified showing a linear and 
\g\ = 2 compatible Zeeman splitting with magnetic field, 
whereas other conductance peaks are difficult to analyze 
similar to investigations in GaAs quantum dots pjJl [20] . 
The observed spin sequence with spin polarization and 
the absence of clear degenerate ground and excited states 
are in agreement with the strength of Coulomb interac- 
tions expected in graphene. 
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bling and T. Heinzel for help and discussions. Support by 
the ETH FIRST Lab, the Swiss National Science Founda- 
tion and NCCR nanoscience are gratefully acknowledged. 



the dielectric constant a mixture of vacuum and Si02 
e r = 2.5 is used. The splitting of the eigenstate energies 
due to the exchange interaction can be estimated using 
£(r s = 1) ~ 0.6A [26-28 . The exchange energy split- 
ting is comparable to the single-particle-level-spacing A 
and can therefore lead to ground-state spin polarization 
S > 1/2 in agreement with our interpretation and obser- 
vations in GaAs quantum dots [T9l [20] . An alternative 
reason for spin-polarization in graphene could arise from 
the valley degeneracy. However, our measurements show 
no indications of four-fold shell filling as it has been seen 
in carbon nanotubes !2T1 l22l [24] . 

Next we describe an attempt to observe the Zee- 
man splitting of excited states in parallel magnetic field. 
Fig. 4(a) shows part of a Coulomb blockade diamond 
measurement at B = T (this is a close-up of Fig. 2(b) in 
Ref . [T6]) . As in Fig. 1(b) several differential conductance 
peaks are visible. Along the white bar at Vb = 12 mV 
the evolution of the peaks is measured in Fig. 4(c). Indi- 
vidual cross sections are plotted in Fig. 4(b) showing the 
change in amplitude and the movement of the peaks. The 
six peaks labeled in Fig. 4(c) are linked to the stronger 
resonances in Fig. 4(a). The magnetic field dependence of 
the two faint lines between resonance 2 and 5 in Fig. 4(a) 
could not be traced in Fig. 4(c). The evolution of the 
peaks can be attributed to two different slopes. For com- 
parison black dashed lines with slopes differing by the 
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